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Abstract. A one-dimensional, time-dependent model is used to study the chemical, 
microphysical, and radiative properties of volcanic aerosols produced by the Mount 
Pinatubo eruption on June 15, 1991. Our model treats gas-phase sulfur photochemistry, 
gas-to-particle conversion of sulfur, and the microphysics of sulfate aerosols and ash 
particles under stratospheric conditions. The dilution and diffusion of the volcanic 
eruption clouds are also accounted for in these conditions. Heteromolecular 
homogeneous and heterogeneous binary H 2 SO 4 /H 2 O nucleation, acid and water 
condensational growth, coagulation, and gravitational sedimentation are treated in 
detail in the model. Simulations suggested that after several weeks, the volcanic cloud 
was composed mainly of sulfuric acid/water droplets produced in situ from the S0 2 
emissions. The large amounts of S0 2 (around 20 Mt) injected into the stratosphere by 
the Pinatubo eruption initiated homogeneous nucleation which generated a high 
concentration of small H 2 SO 4 /H 2 O droplets. These newly formed particles grew 
rapidly by condensation and coagulation in the first few months and then reach their 
stabilized sizes with effective radii in a range between 0.3 and 0.5 ^tm approximately 
one-half year after the eruption. The predicted volcanic cloud parameters reasonably 
agree with measurements in term of the vertical distribution and lifetime of the 
volcanic aerosols, their basic microphysical structures (e.g., size distribution, 
concentration, mass ratio, and surface area) and radiative properties. The persistent 
volcanic aerosols can produce significant anomalies in the radiation field, which have 
important climatic consequences. The large enhancement in aerosol surface area can 
result in measurable global stratospheric ozone depletion. 


Introduction 

Following major volcanic eruptions, the stratospheric 
aerosol layer, known as the Junge layer [Junge et al 1961] 
was greatly enhanced by particles created from the sulfur- 
bearing gas (such as S0 2 ) emissions as well as ash particles 
[Farlow et al., 1981; Turco et al ., 1982; Hofmann and Rosen, 
1983, 1984, 1987; Deshler et al.. 1992; Sheridan et aL , 1992]. 
The significant enhancement of the stratospheric aerosol 
layer can exert an extra forcing on the global climate system 
by altering the radiation balance and can also influence 
heterogeneous chemical reactions and thus global ozone 
abundance. To fully understand the atmospheric effect of 
volcanic eruption, one needs to clarify the chemical and 
physical mechanisms of the formation, evolution, and per- 
sistence of volcanic clouds in the stratosphere. The obser- 
vational and modeling studies of volcanic clouds formed in 
the stratosphere following major eruptions in the past two 
decades thus have been greatly expanded, particularly, due 
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to the recent occurrence of several major eruptions such as 
Mount St. Helens in 1980. El Chichon in 1982. and Mount 
Pinatubo in 1991; however, only a few of model calculations 
to date can compare directly with observations of aerosol 
properties of volcanic clouds in details. This study is to 
develop a volcanic aerosol model based on Turco Toon s 
tracer model to study volcanic aerosols produced by the 
Pinatubo eruption and compare our model results directly 
with observational data collected following the event. 

The eruption of Pinatubo ( 15°N, 120"E) on June 15. 1991, 
resulted in the most massive stratospheric sulfate cloud in 
this century and provided a unique nature laboratory for 
volcanic aerosol studies. The measurements following the 
Mount Pinatubo eruption have supplied large and diverse 
sets of observations. The eruption columns over the vent 
reached in excess of 35 km [Woods and Self . 1992], and the 
densest clouds were initially located between 20 and 30 km 
in tropical regions [McCormick and Veiga . 1992; Deshler et 
al.. 1992; Jdger, 1992]. Total ozone mapping spectrometer 
(TOMS) and stratospheric aerosol and gas experiment 
(SAGE) II satellite sensors monitored the distributions of 
S0 2 amounts and aerosol optical depths, respectively [Bluth 
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et c il . , 1992; McCormick and Veiga , 1992]. The satellite data 
show that the cloud dispersed westward in the easterly phase 
of the quasi-biennial oscillation (QBO) of the stratospheric 
circulation and encircled the Earth within 3 weeks. Its 
meridional dispersion was relatively slow, and the cloud 
remained in a tropical belt extending from about 20°S to 30°N 
for the first few months after the eruption. 

TOMS data also show that about 20 Mt of S0 2 was 
injected into the stratosphere, which is about 3 times the 
amount produced by the 1982 El Chichdn eruption [Bluth et 
al., 1992]. This large emission of S0 2 provided the source for 
sulfuric acid particles in the Pinatubo cloud. Measurements 
of the aged Pinatubo aerosols indicate that 99% were pure 
H 2 S0 4 /H 2 0 solution droplets, confirming the importance of 
homogeneous nucleation of gas-to-particle conversion as a 
volcanic aerosol formation mechanism [Deshler et al 1992; 
Sheridan et al ., 1992]. The observed size distributions of the 
volcanic aerosols in the stratosphere exhibit multiple modes 
[Deshler et al. , 1992; Pueschel et al . , 1992]. Measured modal 
sizes of sulfate particles generated in situ from gaseous S0 2 
by aircraft and balloon samples, and lidar and satellite 
observations are around 0.3—0. 6 ;xm [Pueschel et a /., 1992; 
Deshler et al ., 1992; Russell et al ., 1993b; M. J. Post and 
C. J. Grund, Lidar observations of atmospheric processes 
involving volcanic clouds, preprint, 1993]. Observed strato- 
spheric aerosol extinction from SAGE II multiwavelength 
measurements clearly show the significant difference in size 
and mass between the background and volcanic aerosols 
following the Pinatubo eruption [Thomason, 1992]. Optical 
depths of several tenths {0. 2-0.5) were observed in Pinatubo 
eruption clouds during the first few months [McCormick and 
Veiga , 1992; Stowe et al., 1992; Valero and Pilewskie, 1992; 
Pueschel et al., 1992], which were even larger than those of 
the El Chichon eruption. Such values for the optical depth 
exceed the ambient values by more than l order of magni- 
tude. 

Several detailed model simulations of volcanic clouds 
have been carried out in the last two decades. Turco et al. 
[1983] studied the chemical and physical properties and 
evolution of the stratospheric clouds produced by the erup- 
tion of Mount St. Helens using a one-dimensional (1-D) 
model. They concluded that the general microphysical be- 
havior and structure of the clouds could be adequately 
reproduced with a 1-D model. However, because the sulfur 
chemistry at the time did not involve the fast reaction of 
HS0 3 with 0 2 , a heterogeneous chemical conversion mech- 
anism of S0 2 into sulfate was needed to explain the rapid 
buildup of sulfate mass following the eruption. Pinto et al. 
[1989] used basically the same modd but with updated 
chemistry to examine the effects of stratospheric volcanic 
clouds generated by very large eruptions. They suggested 
that the physical and chemical impacts of volcanic aerosols 
must be self-limiting due to the highly aonlinear growth and 
coagulation processes leading to larger-sized volcanic aero- 
sols with short lifetimes against sedimentation. Turco [1990] 
simulated the El Chichbn eruption cloud with a 1-D model 
and found that homogeneous nucleation was significant in 
the core of the volcanic eruption; aerosol concentrations 
were dominated by newly formed particles. 

Cadle et al. [1976, 1977] and Cadle [1980] simulated the 
global dispersion of several specific vulcanic clouds using a 
two-dimensional (2-D) meridional circiiation model. Capone 
et al. [1983] developed a 2-D model based upon the earlier 


work of Turco et al. [1979a, b] and Toon et al. [1979a, b] to 
study the El Chichon volcanic clouds. These simulations 
reproduced lidar and satellite observations well but inaccu- 
rately depicted the rate of meridional dispersion of the 
clouds. Recently, Tie et al. [1994] simulated El Chichon 
aerosols using a coupled 2-D dynamical/chemical/micro- 
physical model; the resultant aerosol properties are compa- 
rable to that observed following the El Chichon eruption. 
Meanwhile, Bekki and Pyle [1994] also incorporated the 
sulfate microphysics into 2-D chemistry/dynamical model to 
simulate Pinatubo aerosols. They found significant discrep- 
ancies between model simulations and observation in the 
timing of aerosol loading peak and in the magnitude of the 
surface area due to the absence of homogeneous nucleation 
in their model. Some other multi-dimensional modeling 
efforts have been limited to studies of dynamical transport of 
volcanic clouds or the effects of volcanic clouds on atmo- 
spheric circulation and composition [Boville et al.. 1991; 
Brasseur and Granier , 1992; Pitari , 1993; Young et al., 
1994]. 

We use a 1-D model to study the Pinatubo clouds with 
emphasis on understanding the chemical and microphysical 
mechanisms involved in aerosol formation and evolution, 
and on simulations of aerosol microstructure and radiative 
properties. Our model includes gas-phase sulfur photochem- 
istry, gas-to-particle conversion of sulfur, and aerosol micro- 
physics. The dynamical transport of the volcanic eruption 
cloud is parameterized as horizontal dilution and vertical 
diffusion. With updated emission data from Pinatubo erup- 
tions, we use the 1-D model to reproduce the observed 
properties of the volcanic clouds. 

Model Description 

We incorporated the basic chemical and physical pro- 
cesses discussed by Turco et al. [1979a, b] with updated 
chemical reactions, and most of numerical algorithms devel- 
oped by Toon et al. [ 1988], to build up a sophisticated sulfate 
aerosol model for stratospheric conditions. The major com- 
ponents of this model and their interactions through chemi- 
cal and physical processes are illustrated in Figure l. The 
primary source of the undisturbed stratospheric sulfur is 
OCS which gets into the stratosphere mainly through dy- 
namical diffusion processes from the troposphere (Cn/frtvr. 
1976]. It is then decomposed to produce gas phase H : S0 4 
through a number of photochemical reactions (see fable 1). 
During volcanic eruptions, a large amount of sulfur dioxide 
is directly injected into the stratosphere and converted into 
HtS 0 4 by OH oxidation. Then. H 2 S0 4 vapor interacts with 
stratospheric aerosols through nucleation. condensational 
growth, and evaporation. Three types of aerosol are consid- 
ered in the model, namely, pure involatile particles (Aitken 
particles from troposphere or ash grains injected by volcanic 
eruption), pure volatile (H 2 S0 4 /H 2 0) droplets and mixed 
aerosols which have cores of the first type and shells of 
H 2 S0 4 /H 2 0 solution around the cores. The particles of the 
first type are largely transported from the troposphere under 
nonvolcanic conditions but are strongly disturbed by volca- 
nic ash particles. Pure sulfuric acid particles in the second 
type result only from binary homogeneous nucleation. The 
mixed aerosols result from binary heterogeneous nucleation 
of the solid particles of the first type or from coagulation 
between the other two types of aerosols. Depending on the 
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Figure 1 . A schematic diagram of the chemical and micro- 
physical processes treated in our stratospheric sulfate aero- 
sol model. 


supersaturation of H2SO4 vapor, the mixed and the pure 
sulfuric acid particles can expand through condensation or 
shrink due to evaporation. If a mixed aerosol is completely 
evaporated, then its core is released and returned to the first 
type of aerosol. Self coagulation for each type of aerosol and 
heterogeneous coagulation among different types are all 
included in the model. Gravitational sedimentation is also 
treated in the model for each kind of aerosol, serving as a 
major sink for the total aerosol mass. Such sedimentation 
alters the size distribution as larger ash particles fall out. 

Under normal conditions in the stratosphere, homoge- 
neous nucleation may not be a dominant source of new 


particles: however, it becomes extremely important in vol- 
canic clouds because high concentrations of S0 2 can signif- 
icantly increase the production of sulfuric acid vapor and 
thus acid supersaturation. Once substantial homogeneous 
nucleation happens, a large number of very fine pure H 2 S0 4 / 
H 2 0 droplets will be produced within a short time. These 
newly formed particles then grow quickly to form the 
spectrum of the pure acid aerosols or obtain cores by 
coagulating with other two types to transform to the third 
type of aerosols. Our model is specially suitable for studying 
the rapid formation of new particles as in a volcanic plume 
because a new efficient heteromolecular homogeneous and 
heterogeneous nucleation scheme [ Zhao , 1993] is imple- 
mented. 

Throughout this study, a physical height coordinate is 
used for 1-D simulations with 2-km vertical resolution and 30 
layers from ground to 60 km. Implicit numerical scheme and 
family technique with three families (e.g., F { = OCS, F 2 = 
S, SO, S0 2 , and F 3 = S0 3 , HS0 3 , H 2 S0 4 ) are used to solve 
chemical tracers. The vertical profiles of O and H 2 0 con- 
centrations are fixed in accordance with photochemical 
model predictions by Brasseur and Solomon. OH concentra- 
tions are estimated from S0 2 data collected in volcanic 
clouds following the Pinatubo eruption (see next section). 
All the chemical reaction rates are taken from standard 
compilations [DeMore et al., 1991]. Aerosols are different 
from gas tracers because the size distribution must be 
resolved, yielding another physical dimension and additional 
complexity. In order to cover the wide range of aerosol sizes 
with a fair resolution, we specified the size bins by doubling 
particle volume from one category to the next, resulting in 
about 10 volume intervals per decade of size; 40 size bins 
cover the particle radius range from about 0.001 to 10 ^.m in 
the current version of the model. The details of the basic 
formulation and computational algorithms of the model are 
given by Turco et al. [1979a, b]. Toon et a(. [1979a, b]. Toon 
et al. [1988], and Zhao [1993]. 

Model Physical Setting 

To simulate the Mount Pinatubo volcanic cloud, we gen- 
erated a background aerosol field in the tropical region. 


Table 1 . Photochemical Reactions of Sulfur Components in the Stratosphere 



Reaction 

Rate 


R l 

S + O;^ SO + O 

A , = 2.3 x 10 i: 


R 2 

so + 0 . ^ so, + 0 

A, = 2.6 x 10 'V 

Non 1 

} R 3 

SO, + OH + M HSO; + M 

A,}™ = 3.0 x 1 0 - ’ 1 

. n = 3.3* 



A»» = 1.5 x IQ ' 12 

. m = 0 * 

R 4 

HSO, + O, SO, - HO, 

*- 

II 

X 

f 

R 5 

SO 3 + H,0 H,SO, 

A, = 6.0 x Ur 15 


fl 6 

OCS + O SO + CO 

A h = 2.1 x 10 -"e 


P 

OCS + hv ^ S + CO. A < 288 nm 



J 2 

SO, + hv i SO + O. A < 220 nm 




^Three-body reaction with k ^ = (£ 0 [M]/(t + Ao[M]/A x )) x 0.6 {i I M f *1 , where k {) 

*o 300 ( 773 00)-" and k* = k™°\ 77300} . 
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Figure 2. Flow diagram of the volcanic aerosol model. The 
potential effects of volcanic aerosols on global climate and 
stratospheric ozone depletion indicated with dashed line 
arrows are not included in the current model. 


setting the tropopause height at 15 km. We ran the model 
until it reached a steady state and took the ambient aerosol 
spatial concentration distributions in this state as an initial 
value for the volcanic cloud simulations to which the volca- 
nic sources of S0 2 gas and ash particles are added. The 
estimate of total S0 2 injection from the TOMS measure- 
ments is about 20 Mt from the Pinatubo eruption. Observa- 
tions of ash particles and their size distribution in the initial 
volcanic plume are very rare. For model simulations, we put 
20 Mt of S0 2 and 100 Mt of volcanic ash into the model 
stratosphere uniformly distributed between 20 and 30 km. 
Because of absence of observations, the size distribution of 
ash particles used here is assumed similar to that in an earlier 
study by Turco et al. [1983], e.g., a bimodal lognormal size 
distribution for the silicates with 97.5% of the mass in a 
mode at r = 3 ^m with cr = 0.8, and 2.5% of the mass in a 
mode at r = 0.5 jum with cr = 0.5. 

The flow diagram (Figure 2) shows the basic characteris- 
tics of the volcanic aerosol model. After the eruption, the 
large silicate particles quickly fall out off the stratosphere by 
gravitational sedimentation, while the S0 2 from the volcanic 
injection is transformed into H 2 S0 4 by photochemical pro- 
cesses. H 2 S0 4 concentrations can be greatly increased over 
ambient. Once the supersaturation of H 2 S0 4 vapor exceeds 
the critical value for homogeneous nucleation, a large num- 
ber of H 2 S0 4 /H 2 0 solution droplets will be formed by 
heteromolecular homogeneous nucleation. The newly 
formed aerosols can grow rapidly to measurable size be- 
cause of efficient condensation and coagulation and domi- 
nate the volcanic clouds shortly after the eruption. They can 
stay in the stratosphere for several years owing to their 
relatively small size and can travel globally with the atmo- 
spheric circulation, thereby playing a potential role in block- 
ing atmospheric radiation and in catalyzing global ozone 
depletion. 

The chemical oxidation of S0 2 into H 2 S0 4 is crucial for 
aerosol formation and evolution. The conversion rate is 
mainly controlled by OH concentration. Very limited direct 
measurements of OH concentrations have ever been made in 
volcanic clouds. Read et aL [1993] crudely estimated OH 
concentrations as 1.2-3. 3 x 10 6 cm” 3 between 21 and 26 
km in the Pinatubo cloud, using the observed decay rate of 
S0 2 . To reasonably estimate the OH concentrations and 


thus the conversion rates of S0 2 in the Pinatubo cloud, we 
ran the chemistry model for different OH concentrations and 
fit the model outputs with observations. Figure 3a indicates 
the changes in total S0 2 mass with time, associated with 
three different levels of OH (Figure 3b). Comparing the 
model results (curves) with the data (scattering symbols) 
observed by TOMS [Bluth et al 1992], aircraft [ Mankin et 
al. , 1992], solar backscattered ultraviolet (SBUV) [Me- 
Peters , 1992], and MLS [Read et al.. 1993], we found that 
the model prediction with an OH concentration of 2 x 10 6 
cm” 3 below 30 km generally agrees well with the data 
collected from different sources during the first 5 months. 
The predicted total S0 2 mass with this OH level appears 
slightly higher than the MLS measurements after about 150 
days. The inconsistency may result from the fact that S0 2 
mass had fallen below the uncertainty level of MLS mea- 
surement after that time. Thus we use 2 x 10 6 cm 3 of OH 
concentration in the simulation of Pinatubo volcanic aero- 
sols. 

The dynamical processes in this 1-D model are parameter- 
ized as vertical diffusion and horizontal dilution. We use a 
vertical profile of diffusion coefficients similar to that used by 
Turco [1979a, b]. The horizontal dispersion of volcanic cloud 
was simply treated using a cloud area expansion function. 
According to observations [Bluth et al.. 1992; McCormick 
and Veiga , 1992; Stowe et al 1992], Pinatubo volcanic 
plumes initially propagated zonally and encircled the Earth 



Time iday> 


< 30 



OH concentration |10 6 cm' 3 ! 


Figure 3. (a) Changes of total S0 2 masses versus time in 

volcanic clouds; model calculations (curves) associated with 
three different OH profiles as shown in Figure 3b and 
observations (symbols) collected from aircraft and satellite 
following the eruption of Mount Pinatubo. (b) The three OH 
profiles used in calculations for Figure 3a. Read 1.0E + 02 as 
1.0 x 10 2 . 
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within a month. The meridional dispersion is relatively slow 
and it took a year or so to cover the whole globe. For 
simplicity, we assumed that the cloud area expands with 
time as shown in Figure 4, maintaining a uniformly mixed 
composition horizontally. The concentration of each species 
in the cloud at any level is diluted as ambient air is entrained. 
The rates of the dilution of cloud materials are determined by 
the rates of horizontal spreading of the cloud. By considering 
two separate periods, the expansion of cloud area in Figure 
4 is expressed as 

A(t) = A 0 + K x (t - (q), when t 0 < t < t x 

A(f) = A, + K 2 (t - / 1), when t x < t< t 2 
Ait) = A fl , when t > t 2 . 

A denotes the cloud cover area (cm 2 ). K x is the initial 
horizontal diffusion coefficient (cm 2 /s), which can be deter- 
mined from the observed spreading speed of the volcanic 
cloud. Also t is time, and f 0 is the initial time of the model 
(r 0 = 0), corresponding to the initial cross section A 0 = 10 
km 2 of the volcanic plume. A i denotes the cloud area at time 
t x and A a - 51A 0 is the total Earth area. After the eruption, 
the volcanic cloud encircles the Earth in about a month and 
covers about 1/3 of the surface, then disperses poleward and 
covers the globe in about a year. With this kind of cloud 
expansion rate we can easily obtain K\ = 0.5A 0 /day and 
K 2 = 0. 1 A 0 /day by setting t x = 30 days and t 2 = 365 days. 
We have neglected the differences in spreading rates at 
different altitudes actually associated with the dynamics of 
the stratosphere. 

Numerical Simulations 

Using optical particle counter measurements, Hofmann 
and Rosen [1983] studied sulfuric acid droplet formation and 
growth in the stratosphere from the 1982 eruption of El 
Chichon. They suggested two important mechanisms for 
particle formation. One is that sulfate droplets are formed by 
homogeneous nucleation from the gas phase due to en- 
hanced H 2 S0 4 vapor derived from S0 2 injection. The other 
is the rapid growth of the existing particles. Pinto et al. 
[1989] conducted numerical experiments with 1-D aerosol 
microphysical and photochemical models. They confirmed 
the importance of the second view and concluded that 
aerosol microphysical processes of condensation and coag- 
ulation act to produce larger particles and limit the total 
number of particles as the S0 2 injection rate is increased. 
However, the importance of homogeneous nucleation for 
volcanic aerosol formation, as indicated by observations of 
pure sulfate particles in volcanic clouds, has never been 
carefully quantified. Numerical simulations carried out here 
provide a consistent and detailed look at aerosol formation 
and evolution in a well-characterized stratospheric volcanic 
cloud. 

H 2 S0 4 Vapor 

With the physical setting described previously, we inte- 
grate the model for 1000 days. The vertical profiles and time 
series of the simulated H 2 S0 4 are shown in Figure 5. 
H^S0 4 , resulting from the conversion of high initial S0 2 
concentrations in the volcanic cloud, build up very quickly 
to levels sufficient for homogeneous nucleation to occur. For 



Figure 4. Volcanic cloud area expansion adopted in the 
1-D volcanic model. We assume that the initial plume cross 
area is A 0 = 10 7 km 2 . 


instance, H 2 S0 4 mixing ratios exceed 0.05 parts per billion 
by volume (ppbv) between 20 and 30 km during the first few 
days after the eruption, more than 3 orders of magnitude 
higher than ambient H 2 S0 4 levels. About the same order of 
H 2 S0 4 mixing ratio was estimated by Hofmann and Rosen 
[1983] in El Chichdn’s eruption plume. The model simulation 
also shows that such high concentrations of H 2 S0 4 vapor 
cannot be maintained for a long period of time because 
H 2 S0 4 vapor is quickly converted to particles by efficient 
condensation onto preexisting particles, volcanic ash grains, 
and freshly nucleated particles. 

H 2 S0 4 vapor concentrations are determined by produc- 
tion from chemical reactions and loss to particles through 
nucleation and condensation. The time constant for chemical 
production of H 2 S0 4 is equivalent to the S0 2 oxidation 
time, about 30 days, which is much longer than the time 
constant for the loss of H 2 S0 4 vapor to aerosols. Hofmann 
and Rosen [1983], using optical counter measurements at 25 
km, inferred an tv ^-folding" time constant of about 15—20 
min for H 2 S0 4 loss to particles in the El Chichon cloud. The 
preeruption average aerosol distribution at 25 km has a 
surface area sufficient to absorb H 2 S0 4 vapor in about 8 
hours, still substantially faster than S0 2 oxidation. The time 
constant for nucleation is even much smaller than that of 
condensational growth. Therefore although there is contin- 
uous conversion of the volcanic emissions of S0 2 into 
H 2 S0 4 , the acid vapor can not accumulate to high concen- 
trations. Rather, it is almost immediately converted to 
particles. The conversion rates are subsequently controlled 
by the'rates of transformation of S0 2 into H 2 S0 4 (e.g., OH 
concentrations). 

The H 2 S0 4 concentrations in the stratosphere can be even 
lower than their ambient levels about 1 year after the 
eruption. This situation occurs because the production of 
H 2 S0 4 almost returns to the background rate while the 
volcanic aerosol surface area remains enhanced in the low 
stratosphere, creating a greater then normal sink for H 2 S0 4 
vapor. Lower H 2 S0 4 concentrations, in turn, slow particle 
growth. 

Aerosol Size Distribution 

The simulated aerosol number size distributions (Figure 6) 
show that volcanic aerosols are distinct from background 
particles in both concentration and size spectrum. The 
aerosol concentrations are greatly enhanced and the size 
distributions exhibit multiple modes. One of the most signif- 
icant changes is the appearance of very small particles with 
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H 2 S0 4 mixing ratio (ppbv) 



Figure 5. Simulated evolution of H 2 S0 4 mixing ratio 
(ppbv) following the Mount Pinatubo eruption, (a) Vertical 
profiles at different times; (b) time series at different alti- 
tudes. Read IE - 7 as 1 x 10 7 . 

the smallest sizes near 10 _3 The number concentrations 
of these small droplets can exceed 10 5 cm -3 near 21 km. 
Obviously, the high number concentrations of small droplets 
are created in situ through binary homogeneous nucleation 
that is stimulated by high H 2 S0 4 concentration. These 
newly formed aerosols grow rapidly by condensation while 
their number concentrations are reduced significantly by 
coagulation. By 10 days after the eruption, the peak radius in 
number size distribution has been close to 0.1 /xm already, 
with number concentrations of a few hundreds cm -3 . After 
about 1 month following the eruption, new particle formation 
by homogenous nucleation almost ceases when enough 
sulfate mass (and surface area) has been generated to con- 
sume acid vapor effectively by condensation, and thus to 
reduce the supersaturation below its critical value for binary 
homogenous nucleation. Condensational growth and coagu- 
lation then continue to finalize the spectrum of residual 
volcanic sulfate aerosols. After 3 months, less than 50 cm -3 
are left and the condensational growth of these sulfates 
becomes slow due to the reduction of acid vapor concentra- 
tion; aerosol size thus is stabilized at around 0.4 /xm. 

Another distinct feature is the volcanic dust initially 
pumped into the stratosphere with peaks near 0.5 ^m and 3 
/xm. These large ash particles as shown in the figure fall out 
very quickly and most of them disappear after a few months. 
Along with newly formed sulfate droplets, ash grains and 
preexisting ambient particles, which can be heterogeneously 


nucleated by the H 2 S0 4 /H 2 0 vapors, also grow by conden- 
sation. Thus a substantial amount of sulfur vapor is trans- 
formed to particles within the first few hours and days, 
whereas the conversion rate of S0 2 to H 2 S0 4 does not 
increase. Therefore the production of sulfate mass is limited 
by the chemical conversion rate of S0 2 . 

The mass size distributions (Figure 7) also delineate three 
distinct volcanic aerosol modes. The mode at the smallest 
sizes is composed of newly formed sulfuric acid/water drop- 
lets. The other two modes both represent ash particles, 
which initially dominate the total aerosol mass in the strato- 
sphere, but which settle out after about a month. The 
ambient particle mass is too low to be seen on the plots of 25 
and 21 km. The sulfate mode increases in size due to 
condensation and coagulation processes and becomes the 
dominant fraction of the total aerosol mass after ash grains 
fall out. The sulfate mass increases during the first few 
months when condensational growth is significant and then 
decreases when condensational growth becomes slow and 
sedimentation is a dominant sink of aerosol mass. 

The volcanic sulfate aerosols can stay in the stratosphere 
up to several years and significantly alter the nature of the 
Junge layer. The aerosol number and mass size distributions 
after more than 2 years following the eruption show that 
volcanic aerosols are still dominant in the lower stratosphere 
although the background aerosols are recovering. In the 
upper troposphere, the volcanic cloud also significantly 
enhances aerosol mass although it has a minor contribution 
to the aerosol concentrations (see the plots of 15 km). Some 
large-sized particles are added into the troposphere by 
sedimentation from the stratosphere, and changes in mass 
size distribution are obvious. 

The number size distributions, after most ash grains fall 
out, are sharply limited near 1 ^m in the stratosphere. The 
stabilized mode radii in number size distributions are close 
to that in mass size distributions, around 0.4 /xm. This value 
is close to the observed 0.3—0. 5 /xm “effective” radius (e.g., 
the third moment of size spectrum divided by the second 
moment) deduced from advanced very high resolution radi- 
ometer (AVHRR) satellite data and aircraft sampling follow- 
ing the Pinatubo eruption [Russell et al .. 1993a, b; Pueschel 
et al ., 1992]. Figures 8a and 8b show the time series of 
simulated effective radii of total aerosols and pure sulfates, 
respectively. The high initial values of effective radii for total 
aerosols result from large ash particles and decrease sharply 
as the ash grain fall off during the first few weeks. After then, 
the effective radii increase because of the efficient growth of 
new sulfate droplets as shown clearly in Figure 8b. The 
maximum values generally appear around one-half year 
when the condensation growth becomes insignificant. 

Comparing Figures 8a and 8b, we find that the effective 
radii of total aerosols are larger than those of sulfates, 
especially in lower stratosphere during the first 2 years, 
which implies that considerable amount of mass for mixed 
aerosols, sulfates but with ambient or ash particles as cores, 
are still present in the lower stratosphere. However, 2 years 
after the eruption, the situation is opposite. The radii of 
sulfates decrease slower and have larger values than those 
for total aerosols because the larger mixed particles are 
removed by sedimentation, and small Aitken particles are 
transported into the stratosphere by vertical diffusion, 
whereas the sulfates have moderated sizes and do not show 
a significant change in effective radius. 
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Figure 6. Aerosol number size distributions simulated in the Pinatubo eruption cloud with a one- 
dimensional model. Read IE + 5 as 1 x 10 5 . 


From Figure 8b we note that there are different aerosol 
growth rates at different altitudes. For example, the effective 
radius reaches its maximum values around one-half year 
after eruption at 25 km, whereas the maximum appears over 
a year at 19 km. One-dimensional model simulations of 
volcanic clouds carried out by Pinto et al. [1989] indicated 
that the maximum radius for the case with 10 Mt of S0 2 
injection appeared about 1 year with its value of 0.4 ^tm. 
Aircraft measurements also show a similar feature but with a 
slightly larger value of aerosol size, which was probably 
affected by the eruptions of Mount Spurr in Alaska in 1992 
[Goodman et al., 1994]. Hofmann and Rosen [1983] ob- 
served that aerosol growth in the El Chichdn volcanic clouds 
ceased after 7 months. Although the 1-D model results 
cannot be directly compared with observations collected 
from different locations and altitudes, the model simulations 
of the Pinatubo volcanic aerosols are generally reasonable 
and consistent with these earlier modelings and observa- 
tions. 

Vertical Distribution 

The vertical profiles of aerosol concentrations for radii 
greater than 0.15 jim (Figure 9a) show that the Pinatubo 
volcanic eruption greatly enhanced the background aerosols 
in the entire region from the tropopause to about 35 km. The 
densest aerosol cloud is near 20 km with peak concentration 
over 100 cm* 3 within the first few weeks. By about 2 years, 


the aerosol concentration is still significant above its ambient 
value, and the peak altitude is down to 17 km, close to the 
tropopause. The aerosol concentration does not completely 
return to its background level until after 3 years. The vertical 
profiles of aerosol concentration for radii greater than 0.25 
(Figure 9b) show some different features. The peak 
concentration at 10 days is lower than that at 30 days. 
Obviously, most large particles must be ash grains at 10 
days. By nearly 30 days after the eruption, the newly 
formed, initially small sulfates grow to sizes larger than 0.25 
/xm and dominate the aerosol concentration after that. Com- 
pared with Figure 9a, more than 50 percent of the sulfates 
have grown larger than 0.25 jim in radius near 30 days and 80 
percent after one-half year following the eruption. 

Bailoonbome measurements of volcanic aerosols follow- 
ing the Pinatubo eruption at Laramie [De shier et al., 1992] 
show several tens of aerosol concentrations for radii greater 
than 0.15 /x m during July and August. These measurements 
might represent the edge of the volcanic cloud because the 
densest part had not arrived at the observational site, but the 
significant enhancement of aerosol concentration shows 
similar features to our model simulation. The number con- 
centrations of residual volcanic aerosols observed with 
bailoonbome instruments by Deshler et al. [1993] were still 
an order of magnitude higher than their ambient values in the 
lower stratosphere by 1 year after the eruption. The model 
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Figure 7. Aerosol mass size distributions simulated in the Pinatubo eruption plume with a one- 
dimensional model. Read IE — 2 as 1 x 10 2 . 


simulations produce the reasonable magnitude of aerosol 
concentrations and vertical distributions, compared to the 
observed stratospheric aerosol increase following the 
Pinatubo eruption. The enhancement of stratospheric aero- 
sols by the Pinatubo eruption persists for several years. The 
simulated decay rate of volcanic aerosols also indicate 
similar features as those of El Chichdn’s volcanic cloud 
measured by Hofmann et al. [1987]. 

Mass Loading 

The simulated time series of aerosol mass mixing ratios for 
different altitudes is shown in Figure 10. Aerosol mass ratios 
have very high initial values resulting from the ash grain 
input and decrease dramatically as the large ash particles fall 
out. By a few months after the eruption, aerosol mass 
loading in the stratosphere from the volcanic eruption is 
dominated by sulfates (also see Figure 6), and the mass 
mixing ratios still exceed 100 parts per billion by mass 
(ppbm). The observed aerosol mass mixing ratios at 
Laramie, Wyoming (41°N), in Figure 10 indicated that the 
maximum values are above 100 ppbm and appear nearly 180 
days after eruption. After about one-half year, the data fit 
with the model calculations very well although there is 
significant difference between model results and observa- 
tions prior to that time. Obviously, the much lower observed 



Time (day) 

Figure 8. Variations of aerosol effective radii with time 
from 17 to 25 km simulated in the Pinatubo eruption cloud 
with a one-dimensional model: (a) for total aerosols, and (b) 
for pure sulfates (the second type of aerosols). 
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mass mixing ratios (comparing to the model results) within 
the first few months were collected near the edge of the 
cloud because the densest part had not dispersed to middle 
latitude until one-half year. Hence the misfit between the 
model results and observations during that period results 
from inhomogeneities of the volcanic cloud. Nevertheless, 
the model simulation reasonably agrees with the data in 
aerosol mass mixing ratio after the cloud becomes relatively 
uniform. 

Model results also show that the decay rate at higher 
altitude is faster than that at lower latitude. The highest 
values appear near 17-19 km after about 1 year following the 
eruption, which also agrees well with observations from 
Deshler et al. [1993]. Even at 600 days following the eruption 
the magnitude still exceeds 50 ppbm in the lower strato- 
sphere, shown in both simulations and observations. The 
simulated aerosol mixing ratios have still not recovered their 
background values (—0.4 ppbm near 20 km) after 1000 days 
of the eruption. The sulfate mass loading retrieved from 
aircraft measurements shows about an order of magnitude 




Aerosol concentration (cm* 3 ) 


Figure 9. Vertical profiles of volcanic aerosol concentra- 
tions predicted for Pinatubo volcanic eruption cloud (a) for 
particle radii greater than 0.15 ^.m and (b) for radii greater 
than 0.25 fx m. 



Figure 10. Variation in aerosol mass mixing ratios follow- 
ing the Pinatubo eruption, predictions (lines) at several 
altitudes and observations (circles) in the stratosphere 
[Deshler et a/., 1993]. Read IE + 4 as 1 x 10 4 . 

higher than background 2 years after the eruption. The 
model calculations are also consistent with these observa- 
tions. 

Surface Area 

The changes of simulated aerosol surface areas over time 
(Figure 11) show some similar features to the changes in 
mass mixing ratios. The high initial values of simulated 
surface area drop rapidly during the first few weeks because 
of the efficient sedimentation of large ash particles and then 
gradually decay with time. The surface area at 17 km has a 
different behavior during the first few months, because the 
initial cloud was assumed between 20 and 30 km. The 
predictions are generally comparable to the surface areas 
deduced from aerosol sample by Deshler et al. [1993] after 
about 6 months following the eruption. The measurements, 
however, exhibit large fluctuations and lower aerosol surface 
areas within one-half year compared with the model calcu- 
lations. Again inhomogeneities in volcanic clouds and the 
fixed location of the observation site contribute to these 
differences at early times. After 5 to 6 months both measure- 
ments and model results reach a relative constant surface 
area of about 30-60 /xm : /cm 3 between 17 and 21 km, 
although some variations remain in the observational data 
apparently due to transient dynamic transport and seasonal 
processes. In addition, considerable increase of aerosol 
surface area near 500 days appears in the observational data. 
It is probably partially attributed to the dispersion of the 
aerosol loading from eruptions of Mount Spurr in Alaska 
(61°N. 152°W, three eruptions have been reported on June 
27, August 18, and September 16-17. 1992, with injection 
plumes penetrating into the stratosphere (Eos. Transactions , 
American Geophysical Union , volume 74, May 1993)). 

Vertical profiles of the volcanic aerosol surface areas are 
shown in Figures 12a and 1 2b. The volcanic aerosols indicate 
a significant enhancement in surface area from the 
tropopause to 30 km both in simulation and observation. The 
estimation of surface area from lidar and aureolemeter 
measurements in Japan exceed 40 /im : /cm 3 in November 
and December (Figure 12b [Hayashida and Sasano , 1993]). 
The predictions (Figure 12a) of surface areas agree reason- 
ably well with the vertical profiles retrieved from lidar and 
aureolemeter measurements; 40 /xnWcm 3 of surface area 
was also observed about 180 days after the Pinatubo erup- 
tion by balloonborne instrument [Deshler et al 1993]. The 
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Figure 11. Variations in aerosol surface areas following 
the Pinatubo eruption, predictions (lines) at several altitudes 
and observations (circles) in the stratosphere [Deshler et al . , 
1993]. Note the apparent seasonal variations in the observa- 
tions. Read IE + 3 as 1 x 10 3 . 


aerosol surface area increased at least tenfold above back- 
ground values of 1 /um 2 /cm 3 and have remained elevated for 
2-3 years. The enhancement in the aerosol surface area 
possibly causes measurable global ozone depletion in the 
stratosphere. 

Optical Properties 

The significant quantities of stratospheric sulfate aerosols 
generated by Pinatubo substantially perturbed atmospheric 
radiation and temperature. To determine the radiative prop- 
erties of the volcanic aerosols, we use Mie theory to calcu- 
late extinction coefficients and optical depths for the simu- 
lated aerosol size distributions and vertical profiles as a 
function of time. The refractive index of sulfate aerosols is 
temperature and concentration dependent, but we use a 
constant value of 1 .45-0/. The choice of refractive index has 
little effect on the calculated optical depth for the range of 
values that might be expected assuming either sulfuric acid 
or dust particles [Russell and HamilU 1984]. 

The observational data show that the optical depth is 
wavelength dependent during the first month and becomes 
relative flat after 2 months (Figure 13a). The model optical 
depth (Figure 1 3b) also show distinct wavelength depen- 
dence within the weeks after the eruption, especially, on the 
day-one curve. After a few months, the model optical depth 
becomes much flatter (with less wavelength dependence). A 
significant wavelength dependence of the optical depth indi- 
cates the presence of small particles in the volcanic cloud. 
This relationship agrees with the evolution of the aerosol 
size distributions in the volcanic cloud discussed in the 
previous section. The model optical depths in the early time 
in Figure 13b exceed the observed values in Figure 13a 
because the early observations at Mauna Loa were sampling 
the edges of the dispersing volcanic clouds. 

The calculated optical depths of the volcanic aerosols at a 
wavelength of 0.5 ^im (Figure 14) indicate an increase of 1 to 
2 orders of magnitude from ambient values (—0.005) in the 
stratosphere. The simulated cumulative optical depths are in 
the range of 0.4— 0.5 during the first few months, resulting 
initially from ash and then from sulfates after about a few 
weeks. The vertical profile at 90 days shows the highest 
values among the others. The major contribution to this large 
value of optical depth is mainly contributed by sulfate 
particles because most ash has been removed from the 
stratosphere by that time. By a year after eruption, optical 


depths of volcanic aerosols are still as high as 0.1 (see the 
profiles at 360 days), and near 0.05 at 600 days at altitudes 
near 20 km. Even by 1000 days after the eruption, the optical 
depths have not completely returned to their background 
values. 

The time series of cumulative optical depth at different 
altitudes in the stratosphere is plotted in Figure 15. During 
the first few weeks, the optical depths appear up and down 
twice, resulting from the complexity of multimode size 
distributions including both ash particles and small sulfates 
and of spectrum changes due to removal and growth of 
different sized aerosols. By a month after the eruption, most 
contributions on optical depth are from sulfates. The cumu- 
lative optical depths reach their maximum values around 2-3 
months and the values range between 0.14 and 0.55 from 25 
km to 17 km, and then decay with time in a consistent way 
as aerosol mass and surface area. 

There are numerous measurements about volcanic aerosol 
optical properties observed at different wavelengths, alti- 
tudes, and latitudes following the Pinatubo eruption. The 
stratospheric optical depth in the months following the event 
was observed in the range between 0.2 and 0.5 [Stowe et al ., 
1992; Valero and Pilewskie, 1992; Russell et al ., 1993a, b]. 




Figure 12. (a) Simulated vertical profiles of aerosol surface 

areas at different times after the Mount Pinatubo eruption, 
(b) The aerosol surface area of stratospheric aerosol on 
November 25, December 1, and December 13, 1991, esti- 
mated from the lidar and aureole meter measurements at 
36.0°N, 140. 1°E [Hay as hid a and Sasano , 1993]. 
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The optical thickness derived from reflected solar radiation 
measurements (at A = 0.5 fim) by the advanced very high 
resolution radiometer (AVHRR) on the NOAA 1 1 satellite 
show peak optical depths exceeding 0.5 at some locations 
during July and August 1991. A mean optical thickness of 
0.31 was estimated on August 23 [Stowe et al ., 1992]. 
Midvisible optical depths higher than 0.4 were observed 
from July 7 to 14, 1991 [Valero and Pilewskie, 1992]. The 
optical depth measured at Mauna Loa observatory by an 
Ames Research Center tracking sunphotometer was 0.2 
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Figure 13. Total aerosol optical depths versus wavelength, 
(a) Observations at Mauna Loa observatory at different 
times following the Mount Pinatubo eruption on June 15, 
1991. The June 1991 data roughly correspond to background 
conditions but include substantial tropospheric aerosols. 
The July data show a strong wavelength dependence. The 
optical depths had a flat spectrum by August/September and 
were reduced about half value by July 1992 [Russell et al . , 
1993b]. (b) Model results, giving the optical depth spectrum 
for the total optical depth above 21 km. Unlike the data tn 
Figure 13a, tropospheric aerosol effects are not included. 
NASA American Association for the Advancement of Sci- 
ence, solid line; National Oceanic and Atmospheric Admin- 
istration CMDL, dashed line. Read IE + 0 as 1 x 10 . 



Accumulative aerosol optical depth 

Figure 14. Vertical profiles of optical depth for the simu- 
lated Pinatubo volcanic clouds at different times after the 
eruption. The ambient optical depths are also shown for the 
background aerosol model. The aerosol optical depths are 
cumulative, from the top of the atmosphere. Hence the 
cumulative optical depth at each level represents the total 
optical depth of the aerosol lying overhead. 


about 2 months after the Pinatubo eruption [Pueschel et al.. 
1992; Russell et al., 1993a]. In March 1992. optical depths as 
large as 0.22 were still observed over large areas in the 
stratosphere [Russell et al., 1993a], The global average 
optical depth of Pinatubo aerosols (after the cloud dispersed 
over to the whole globe) was estimated about 0.15 in early 
1992 [Satoetal., 1993 ]. Although our 1-D model calculations 
can not precisely compare with these measurements from 
different locations and altitudes collected with a variety of 
instruments, the basic features and magnitude of aerosol 
optical depth between model results and observations are 
very close and consistent. The significant increase in the 
atmospheric optical depth can produce measurable pertur- 
bations to global temperatures and even the general circula- 
tion of the stratosphere and troposphere. Stratospheric 
warming and surface cooling were observed during the year 
following Pinatubo [Angetl. 1993]. Hansen et al. [1992] and 
Graf et al. [1992] predicted using general circulation model 
experiments that significant climatic impacts would result. 



Time (day) 

Figure 15. History of cumulative optical depths of 
Pinatubo volcanic aerosols predicted with the present model 
at different altitudes. 
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and these have materialized in the form of a peak global 
cooling of about 0.5°C. 

Conclusions 

A one-dimensional volcanic cloud model, including chem- 
ical and physical processes involving sulfate aerosols, was 
developed to study the formation and evolution of the 
Pinatubo eruption clouds in the stratosphere. By comparing 
model simulations to observations of the Pinatubo volcanic 
cloud, the model has been shown to reasonably reproduce 
the basic aerosol properties, including their size distribution, 
total mass and concentration, temporal and altitude varia- 
tions, and radiative properties. The model, however, cannot 
reproduce the latitudinal variations observed throughout the 
several year evolution of the Pinatubo cloud. 

The model includes all of the important aerosol formation 
processes in volcanic clouds. The S0 2 injected by the 
volcanic eruption is transformed into sulfuric acid vapor by 
OH oxidation. The enhanced H 2 S0 4 concentration creates 
new H 2 S0 4 /H 2 0 particles through homogeneous nucle- 
ation. These fine droplets then evolve through condensa- 
tional growth and coagulation. Model results show that the 
homogeneous nucleation of H 2 S0 4 /H 2 0 droplets occurs 
within the first few weeks. Aerosol concentrations are sub- 
sequently controlled by coagulation. Condensational growth 
of the sulfate particles is altitude dependent and generally 
becomes insignificant about one-half year after the eruption. 
The effective radius of the stabilized volcanic sulfate aero- 
sols is around 0.4 p,m. 

Clearly, the late-time density of a stratospheric eruption 
cloud is mainly determined by the sulfur gas emission, rather 
than the ash emission. The ash particles, which are quite 
large, fall out very quickly. The residual volcanic clouds are 
then composed mainly of H 2 S0 4 /H 2 0 droplets generated in 
situ from the sulfur emissions. Homogeneous nucleation is 
essential for the production of new sulfuric acid particles, 
which are observed in the stratosphere. Our model calcula- 
tions agree with observations indicating that more than 90% 
of the volcanic aerosols are completely volatile. It is the 
sulfate particles that reside in the stratosphere for several 
years and spread over the planet, greatly intensifying the 
Junge layer. The 2-D model study by Bekki and Pyle [1994] 
also indicated the importance of homogeneous nucleation. 
Without including homogeneous nucleation, the simulated 
timing of aerosol loading peak and the magnitude of the 
surface area exhibit significant differences from observa- 
tions. 

Simulated aerosol surface areas and mass mixing ratios 
within 2-3 years of the eruption are also comparable to those 
observed at 4l°N by Deshler et al. [1993], but the observa- 
tions indicate important variations associated with the lati- 
tudinal inhomogeneity of real volcanic clouds. Many of the 
measured parameters have lower values than predicted by 
the model during the first few months. Once the clouds have 
become relative uniform in the stratosphere, the agreement 
of the model with point observations improves. The aerosol 
surface areas in the lower stratosphere are still more than 10 
/xm 2 /cm 3 , comparable to the surface areas of polar strato- 
spheric clouds, even 2-3 years after the Pinatubo eruption. 
This large enhancement in the aerosol surface has apparently 
led to measurable global ozone depletion in the stratosphere 
[Grant et al., 1992; Hofmann et al., 1993; Schoeberl et al . , 
1993]. 


Corresponding to the evolution of the aerosol size distri- 
bution in the volcanic clouds, the optical depths show a 
wavelength dependence within the first lew weeks associ- 
ated with a large number of small particles generated by 
homogeneous nucleation. A few months after the eruption, 
the wavelength dependence of the optical depth has become 
relatively flat, because the sulfate particles have grown to 
larger sizes through condensation and coagulation. In gen- 
eral, optical depths associated with volcanic sulfate aerosols 
in the model simulations remain at least an order of magni- 
tude above ambient values during the first 2 years following 
the sulfur injection. Our model predictions are consistent 
with a variety of optical measurements following the erup- 
tion of Mount Pinatubo. The significantly enhanced optical 
depth would be expected to cause, and apparently has 
resulted in, observable global temperature perturbations as 
well as major chemical changes throughout the stratosphere. 

The observed variability in aerosol properties can only be 
simulated with a multidimensional model, preferably in three 
dimensions. Such a model would allow the full range of 
observable conditions to be sampled numerically. These 
conditions include spatial structures (e.g.. vertical and hor- 
izontal), microphysical structures (e.g.. size distribution and 
optical depth variations), and chemical structures (e.g.. SO : 
distribution and decay). Nevertheless, the genera! character- 
istics and long-term variations of the key parameters seem to 
be well reproduced by the present one-dimensional model, 
when an appropriate initialization is possible. Given that the 
present model treats all of the physical and chemical pro- 
cesses thought to contribute to the evolution ot aerosol 
properties, although not the dynamical processes, and can 
generate all of the ‘'observables available trom measure- 
ment systems, the next logical step, which we plan to take, 
would be to couple the volcanic aerosol model to a three- 
dimensional model. 
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